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Abstract 

The coupling constants g pm and g ur)1 are calculated using QCD sum rules 
method by studying the three point prj'j and W777 correlation functions. A 
comparison of the results with the values of the coupling constants that are 
deduced from the experimentally measured decay widths of p — > 7/7 and uj — > 
777 decays is performed. 
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The method of QCD sum rules is one of the efficients tools for studying hadron physics. 
This method has been successfully applied to calculate many hadronic observables, such as 
decay constants and form factors [1,2]. On the other hand, radiative transitions between 
pseudoscalar (P) mesons have been an important area of study in low-energy hadron physics 
for more than three decades. These transitions have been analyzed within the frameworks 
of phenomenological quark models, potential models, bag models, and also by employing 
effective Lagrangian methods [4,5]. The radiative transitions V — > P7 are characterized by 
the coupling constants gv rn - Since low energy hadron physics is governed by nonperturbative 
QCD, it is very difficult to obtain the numerical values of these coupling constants from the 
first principles. For this reason, some specific nonperturbative methods have to be developed 
to be used as calculational tools. Among these methods QCD sum rules have proved to be 
very useful to extract the coupling constants. A recent review of QCD sum rules method is 
provided in [6] where more references can also be found. 

In this work, we calculate the coupling constants g PV7 and g ur)1 associated with the 
radiative decays p — > 777 and u — * 777 by employing the traditional QCD sum rules method 
which provides a model independent way to calculate the coupling constants. The coupling 
constant g pvl was previously calculated by T. M. Aliev et al. [7] in the framework of light 
cone QCD sum rules. Our analysis, therefore, complements the results obtained in that 
paper. 

In accordance with the general strategy of QCD sum rules method, we begin by consid- 
ering the three point correlation function 

IWp,p') = / WV^ 1 < o\T{fi(0)fi(x)jr,(y)}\0 > (1) 

where the interpolating currents for vector meson p and u are = -^{u^ v u — dj u d), 
j% = -^{u^ v u + drj v d\ respectively. We take 77 — 77' mixing into account and use the 
interpolating current for 77 meson as j v = ■^=(ui f y 5 u + di r y 5 d)cos9 — (5775s) sinO where 9 is 
the mixing angle in the quark-flavour basis. The electromagnetic quark current is given as 
j 1 = e^a^fjU + edd'y^d, where e u and denote the quark charges. 
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The theoretical part of the sum rule for the coupling constant gy vl is calculated by 
considering the perturbative contribution and the power corrections from operators of dif- 
ferent dimensions to the three point correlation function. In the spirit of QCD sum rules 
techniques, we consider the three point correlation function in the Euclidian region defined 
by p 2 = —Q 2 ~ —1 GeV 2 , p' 2 = —Q' 2 ~ —1 GeV 2 . In this region, the perturbative 
contribution can be approximated by the lowest order quark loop diagram shown in Fig. 1. 
Moreover, we consider the power corrections from operators of different dimensions, propor- 
tional to terms <qq>,< qa-Gq > and < (qq) 2 >. Since the gluon condensate contribution 
proportional to < G 2 > is estimated to be negligible for light quark systems, it is not taken 
into account. We perform the calculations of the power corrections in the fixed point gauge 
[8]. Moreover, we work in the SU(2) flavour context with m u =m ( i=m q and we work in the 
limit m q = 0. In this limit, the perturbative quark-loop diagram does not make any contri- 
bution, and only contributions result from the operators of dimensions d=3 and d=5 that 
are proportional to < qq > and < qa ■ Gq >, respectively. The relevant Feynman diagrams 
for the calculation of the power corrections are shown in Fig. 2 and Fig. 3. 

We then calculate the three point correlation function U^ u (p,p') using phenomenological 
considerations. This function satisfies a double dispersion relation. We choose the vector 
and pseudoscalar channels and by saturating this dispersion relation by the lowest lying 
meson states in these channels the physical part of the sum rule is obtained as 

where the contributions from the higher states and the continuum are shown by dots. The 
overlap amplitudes for vector and pseudoscalar mesons are < 1 | V" >= XyUy where Uy 
is the polarization vector of the vector meson V = p, u> and < 77 1 1 >= A,,. The matrix 
element of the electromagnetic current is given by 

< V(p)\j;\ V (p f ) >= - t ^g Vrjl K(q 2 )e^p u u a qp (3) 

my 

where q = p — p' and K(q 2 ) is a form factor with K(0) = 1. This matrix element defines the 
coupling constant gy m through the effective Lagrangian 



C = —9vm^ a %V v d a Apri (4) 
my 

describing the Vr}^- vertex [9]. 

After performing the double Borel transform with respect to the variables Q 2 and Q' 2 , 
we obtain the sum rule for the coupling constant gy r?7 in the form 

9vrn = ——e^e^ 1 [e u <uu> ±e d < dd >) 

( 3 5 2 1 3 2 1 \ 

1 rn n —— m n T cos9 (5) 

V 2 16 °M 2 16 °M' 2 ) V ; 

where the relation < qa ■ Gq >= m\ < qq > is used. In this expression the plus sign is 
for p meson and the minus sign is for u meson. In the numerical evaluation of the sum 
rule the values m 2 = (0.8 ± 0.02) GeV 2 , < uu >=< dd >= (-0.014 ± 0.002) GeV 3 
[6], and m p = 0.77 GeV, m w = 0.781 GeV, m T] = 0.547 GeV are used [10]. The 
overlap amplitudes for vector meson states are calculated using the experimental leptonic 
decay widths of V — > e + e~ decays [10] and the values X p = (0.17 ± 0.03) GeV 2 and 
A w = (0.15 ± 0.02) GeV 2 are obtained. The overlap amplitude for rj meson state was 
determined earlier by QCD sum rules analysis as = (0.23 ±0.03) GeV 2 [11]. We use the 
value of the mixing angle as 9 = —19° ± 2° [11]. 

The dependence of the coupling constants gv V j on the Borel parameters M 2 and M' 2 are 
analyzed by studying the independent variations of M 2 and M' 2 in the interval 0.6 GeV 2 < 
M 2 ,M' 2 < 1.4 GeV 2 since these limits determine the allowed interval for the vector channel 
[12]. We show the variation of the coupling constant g pm and g wrry as a function of the Borel 
parameters M 2 for different values of M' 2 in Fig. 4 and in Fig. 5, respectively. These 
figures indicate that the sum rule is quite stable with these reasonable variations of M 2 and 
M' 2 . We choose the middle value M 2 = 1 GeV 2 for the Borel parameter in its interval of 
variation and obtain the coupling constants gvrn as IWy — 1-2 ± 0.3 and g W7r7 = 0.4 ± 0.06 
where the uncertainties result from the variations of M 2 and M' 2 and from the estimated 
values of the vacuum condensates. 

If we use the effective Lagrangian given in Eq. 4, then the decay width for V — > 777 is 



obtained as 



a {m 2 v - m 2 v f 



r(V — > — — - ~^^~9vin <«) 

We then utilize the measured decay widths T(p — > 777) = (57 ± 10) keV and T(uj — > 777) = 
(5.5 ± 0.9) keV [9] and obtain the coupling constants gvq-y as g pm = 1.42 ± 0.12 and g W7r7 = 
0.42 ± 0.03. Our results are, therefore, in good agreement with the coupling constants 
deduced from the experimental values of the respective decay widths. Moreover, our result 
for g pm is also consistent with the value g p ^ 7 = 1.42 ±0.2 calculated by T. M. Aliev et al. [7] 
in the framework of light cone sum rules, thus our study employing traditional QCD sum 
rule method supplements the previous light cone QCD sum rules calculation. 
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FIGURES 
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FIG. 3. Operators of dimension 5 corrections proportional to < qa ■ Gq >. The dot lines denote 
gluons. 
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FIG. 4. The coupling constant g p „ 7 as a function of the Borel parameter M 2 for different values 
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FIG. 5. The coupling constant g wr?7 as a function of the Borel parameter M 2 for different values 
of M' 2 . 
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